Vortices and chirality in multi-band superconductors 
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We investigate some significant properties of multi-band superconductors. They are time-reversal 
symmetry breaking, chirality and fractional quantum flux vortices in three-band superconductors. 
The BCS (Bardeen-Cooper-Schrieffer) gap equation has a solution with time-reversal symmetry 
breaking in some cases. We derive the Ginzburg-Landau free energy from the BCS microscopic 
theory. The frustrating pairing interaction among Fermi surfaces leads to a state with broken time- 
reversal symmetry, that is, a chiral solution. The Ginzburg-Landau equation for three-component 
superconductors leads to a double sine-Gordon model. A kink solution exists to this equation as 
in the conventional sine-Gordon model. In the chiral region of the double sine-Gordon model, an 
inequality of Bogomol'nyi type holds, and fractional- tt kink solutions exist with the topological 
charge Q. This yields multi-vortex bound states in three-band superconductors. 



I. INTRODUCTION 

Since the discovery of oxypnictides LaFeAsOi_a;Fa;[l[, 
BaFeaAsaia], LiFeAs[l,|i] and Fei+:rSe[i,@, the Fe pnic- 
tides high-temperature superconductors have attracted 
extensive attention. There are numerous experimental 
studies regarding the electronic states of the new fam- 
ily of iron-based superconductor 0-113 ■ The undoped 
samples exhibit the antiferromagnetic transition!^ llol |. 
and show the superconducting transition with electron 
doping The band structure calculations indicate that 
the Fermi surfaces are composed of two hole-like cylin- 
ders around F, a three-dimensional Fermi surface, and 
two electron- like cylinder around M for LaFeAsOflS]. 
This family of iron pnictides is characterized by multi 
Fermi surfaces, and theoretical studies have been based 
on multi-band models with electronic interactions [13- 
[l6j . An importance of multi-band structure is obvi- 
ously exhibited in recent measurements of the Fe isotope 
effect [TtI, [Tsj ■ The inverse isotope effect in (Ba,K)Fe2As2 
can be understood by the multi-band model with com- 
peting inter-band interactions [l^. The two-gap theory of 
superconductivity has a long history, and is the general- 
ization of the BCS theory to the case with two conduction 
bands HS-il 

The objective of this paper is to study time-reversal 
symmetry breaking and fractional quantum-flux vortices 
in multi-component superconductors. We investigate 
physical properties of the superconducting state with 
time-reversal symmetry breaking. The 3x3 gap equa- 
tion has such a solution if the pairing interactions sat- 
isfy some conditions. To investigate vortex solutions, we 
derive the Ginzburg-Landau functional (23l - l25j from the 
microscopic theory of multi-component superconductiv- 
ity. The importance of phase dynamics has been pointed 
out previously [m, I26l430j . The phase variables of the 
order parameters may lead to a new state as a mini- 



mum of the Ginzburg-Landau potential. If the phase of 
the gap function takes a fractional value other than or 
TT, a significant state, called the chiral state[3l|, H^l, ap- 
pears and the time reversal symmetry is broken. This un- 
conventional state is induced from the Josephson terms 
with frustrating pairing interactions. The signs of the 
Josephson terms play an important role in determining 
the ground state. 

The three-band model leads to the double sine-Gordon 
model as a model to describe the dynamics of phase vari- 
ables. This model is regarded as a generalization of the 
usual sine-Gordon model[33j and also has kink solutions 
as in the sine-Gordon model. A new feature that ap- 
pears first in the double sine-Gordon model is that we 
have fractional- TT kink solutions in the chiral. We can de- 
fine the topological charge from the topological conserved 
current. In the chiral region, fractional flux vortices may 
exist on a domain wall of the kink. 

The paper is organized as follows. In Section II, we 
present the model Haniiltonian that is considered in 
this paper. In Section III, we show the model that ex- 
hibits time-reversal symmetry breaking. In Section IV, 
the Ginzburg-Landau functional for three-band super- 
conductors is derived by using the Gor'kov method. Each 
term of this functional is expressed in terms of the matrix 
G = (gij ) where gij is the pairing interaction between the 
bands i and j. Here we investigate the magnetic prop- 
erties near the upper critical field. In Section V we con- 
sider the phase dynamics of the order parameters. We 
investigate the ground state of the potential composed 
of phase variables and show the existence of the chiral 
region with time reversal symmetry breaking. The dou- 
ble sine-Gordon model appears in the three-band model 
and exhibits unique properties that are not contained 
in the conventional sine-Gordon model. In Section VI 
we discuss the existence of fractional flux vortices and 
multi-vortex bound states. We also discuss applications 
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to multi-band superconductors such as Fe pnictides. 

II. THREE-BAND BCS MODEL 

We adopt the three-band BCS model with the attrac- 
tive interactions 

H = J2 [ dv4,{v)K,{v)^U^) 

(1) 

where i and j (=1,2,3) are band indices, ifi(r) stands 
for the kinetic operator. We assume that gij = g*^. The 
second term is the pairing interaction and gij are coupling 
constants. The mean-field Hamiltonian is 

Hmf - E / rfr[E^l(i-)if.(r)V'..(r) 

+ A,(r)4(r)v4^(r) + A*(r)^a(r)V'.t(i-)] , (2) 
where the gap function in each band is defined by 



and its complex conjugate is 
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We define Green's functions as follows [33|. 

Gjaa'{x - x') = -{TTi^ja{x)lp]^,{x')) , 



(3) 



(4) 



(5) 



F+^, (x - x') = {Tr^Ux)^l, (x')), (6) 



where is the time-ordering operator and we use the 
notation x = (r, r). In terms of the Green's functions, 
the gap functions satisfy the system of equations 



A:(r) = E-9*.^itt(^'=^ + 0;^'^) 



This yields the gap equation. 
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2T J ' 



(7) 



(8) 



where Ej = y + l^jP ^-^d T is the temperature where 
we set Boltzmann constant fc^ to unity. Nj is the den- 
sity of states at the Fermi surface. Since all the bands 



couple with each other through mutual interactions gij, 
we have one critical temperature TcpO]. At the critical 
temperature T — Tc, this equation reads[23] 



A,: In 



2e'^a;, 



(9) 



for the cutoff energy Wc- 7 denotes the Euler constant. 
Here we assume the same cutoff energy in all the inter- 
actions. The critical temperature is obtained as a 
solution to the equation 



= 



rigiiNi -~ 1 TOi2A^2 '7513^3 
V921N1 77322^^2 - 1 V923N3 
VdsiNi rigz2N2 VdsaNa - 1 



(10) 



where ry = hi{2e'^uJc/T^Tc). The system of equations in 
eq-© yields a set of differential equations (23il 



A*(r) 



In 



ttT 

7C(3) 
48(7rT, 

7C(3) 



2 E dji^ive 



8(7rr, 



2 Z-^ 



g*,mA*,ir)\Aeir)\' 



A,*(r) 
(11) 



Here, e is the charge of the electron, and vg is the elec- 
tron velocity at the Fermi surface in the £-th band. The 
Planck constant h has been dropped except in front of 
the vector potential A. We set amn = gmnNn- Then, 
the equation for Aj is 







flj-jln 



ttT 



A,+ln( }2 %^A,-t-. 

(12) 

In the second term of the right-hand side we can replace 
T by Tc near the transition temperature. From eq. © we 
obtain 



E 




bii bi2 bi3 

^21 ^22 ^23 
^31 ^32 ^33 




(13) 



where bji are elements of the inverse of the matrix A 
{cLij)'- = (bji). Then the equation for Aj reads 







ajjln 
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A, 
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For example, for j = 1 we obtain 

'2e''ujr\ 1 



(14) 
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(G-i)i3A5 



(15) 



where G = (5^ ) is the matrix of coupling constants. To 
obtain the multi-band Ginzburg-Landau functional, we 
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multiply eq.([T4]) by A*iVj and take a summation with re- 
spect to j. We use the gap equation A| = g^-NjA* 
at T = Tc- Then the energy functional density / is 



9n 



(gN) 



iV,ln^-(G-i) 



ttT 



J3 



+ E|^a;(g-),,a. 

^ ^ 

, 7C(3) 



167r2T2 



(16) 



Here we set 77 = 1/(17 TV), and we neglect unimportant 
constants gjjNj/ (gN) in the following. The fourth order 
term is simply given by (|A^p)^. Please note that out 
functional is relevant near Tr- 



ill. TIME-REVERSAL SYMMETRY 
BREAKING 

The three-band superconductor exhibits a significant 
state with broken time-reversal symmetry, that is not 
shown in single-band and two-band superconductors. 
This originates from the fact that the gap equation for the 
three-band model has complex eigen-functions in some 
cases. The gap equation is 



(giiNi gi2N2 gisNs 
921N1 922^2 523^3 
531^^1 532^^2 533^3 




, (17) 



where — \\i{2e^ijJc/'!TTc). We assume that the matrix 
(gijNj) is real and symmetric. When two eigenvalues 
are degenerate, we obtain two complex eigenvectors by 
making linear combinations. For example, let us consider 
a matrix A = {Aij) given by 



All 


= J(cos^ 


6 cos^ (fi + sin^ (f) + D 


(18) 


A12 


= ^21 = 


— J sin^ 6 sin ip cos (fi 


(19) 


Ai3 


= ^31 = 


— J sin 6 cos cos ip 


(20) 


A22 


= J(cos^ 


9 sin^ ip + cos^ ip) + D 


(21) 


A23 


= ^32 = 


— J sin 6 cos 9 sin p 


(22) 




= J sin^ (. 




(23) 








(24) 



where J is a positive constant and D is also a constant 
satisfying J + D > 0. J indicates the strength of inter- 
band coupling constants. The degenerate eigenvalue is 
given by 



X = J + D. 



(25) 



9 and p are Euler angles in the range oi < 9 < n and 
0<ip <2'K. 



The eigenvectors of the matrix in ea.p4)) that corre- 
spond to the eigenvalue X = J + D are 



(26) 



/ COS 9 cos p> \ 




' sin p \ 


cos 9 sin p} 




— cos p 


\ -sin 9 i 




{ i 



When p> — 7r/4, they are 

— i(p 



1 



1 



(27) 



where the angle (f) satisfies cos (j) = — 1/(\/2 tan6'), and we 
assume that |-\/2tan0| > 1. In this way the time- reversal 
symmetry broken state is obtained from the gap equa- 
tion. If we set p) — Txj^ and tan6' = \/2, we obtain the 
case with three equivalent interband pairing interactions, 
that is, the three equivalent bands[35]. In this case we 
have = 27r/3, and thus the phase difference between 
order parameters of three bands is 27r/3. 



IV. THREE-GAP GINZBURG-LANDAU 
FUNCTIONAL 

In this section, we derive the Ginzburg-Landau free 
energy from the BCS model to investigate chirality and 
kinks. From the results in the section IE, the Ginzburg- 
Landau functional for three-gap superconductors has the 
form 



F 



[dr ^a,|V,f + 

U=i j=i 

(7l2V'tV'2 +72lV^2V'l) 

(723V'2'03 + I32tp3lp2) - (731-03 V"! + 7l3V't'03) 
3 



E^. 



V-f I— A 

00 



0. 



1 



/3j are constants and 0o is the fiux quantum, 



he 



he 
2H 



The coefficients of bilinear terms are 



iV.ln 



ttT 



33 



li3 =~iG 



(28) 



(29) 



(30) 



(31) 



for z, J = 1, 2, • • • . The formula for more than three bands 
is straightforwardly obtained in terms of G~^. The crit- 
ical temperature is given by, with 77 in eq. (|10p . 



(32) 



4 



For the two-band case, if we assume that gap functions 
are constant, 77 is exphcitly written as 



To have a nontrivial solution (Ci, C2, C3) 7^ 0, the secular 
equation must hold. 



1 



ail + 022 - V (an ^ ^22)^ + 4ai2a2i 



2detA 

gii/iV2 + g22/iVi 
2detG' 



ai + Ki/e 
-111 
-731 



—712 —713 
^2 + ^2/^^ -723 
-732 as + i^s/C^ 



0. (43) 



1 /I /gii g22 
detC?V4V^2 A^i 



912921 

N1N2' 



(33) 



Tc and the upper critical field Hc2 — 0o/(27r^^) is de- 
termined from this third-order secular equation. In the 
two-band case this equation is simply 



This is the formula first obtained by Suhl et al. [2(j| In the 
case where three bands are equivalent, we can assume 
simply that gu = g, Ni ^ N{0) and gij = ly (i,j=l,2,3). 
In this case we have 



ai/Ki + 1/^2 -112/ Ki 
-721/X2 012/ K2 + 1/^2 



= 0. 



(44) 



a, 



g + v 



{g + 2ri){g~u) 



(34) 



(35) 



The Hc2 and the coefficients (Ci,C2) are immediately 
derived from this as 



ig + 2iy)ig-iy)' 

Here we investigate the magnetic properties near the 
upper critical field Hc2- The Ginzburg-Landau equations 
are 

aitpi-'-fi2i^2-Ji3'ip3+l3i\ipi\'^i'i-Ki ( V + i^A ) ifji = 0, 



1 / Ql ^2 

2 l^^i X2 




Q!2 



712721 



2i4:i 2K2 J K1K2 

(45) 



Ci 



El 

712 




0^2 



V2ii:i 2K2 



712721 



(36) 



a2^/'2-721'01-723V'3+/52|'02| '02--f^2 V + i — A -02 = 0, 



27r 



00 



27r 



i^ii^2 
(46)" 

Let us consider a simple case for the three-band model, 
where we assume that three bands are equivalent: a.j — 
a, 7„ = J(> 0), Kj = K, gjj = g, g,j = r]{> 0) {i ^ j) 
and Nj — N{0). We easily obtain 



(37) 



Hc2 — 



a3'03-731'01-732V'2+/33|'03| '03--?^3 + ^'^^j "^3 ^ 

(38) 

Let us assume that the magnetic field H is along the 
z axis: H = {0,0, H). We adopt the vector potential 
A — {0,Hx,0). If we suppose that dipj/dy — 0, we 
obtain for j=l,2,3. 



2ttK 



(-a+2 J) = 



2ttK 



iV(0)ln 



2eTwc 
ttT 



1 



Jc = exp 



1 



5 + 277 J 
(47) 

(48) 



li^j- E 7,^'/'^-^m+^.^^f'^'V', =0. (39) 



7V(0)(g + 2,y) 

In the n-band model, these formulae are generalized to 

00 



2ttK' 



-a+{n- 1)J], 



Since this equation is analogous to that for the harmonic 
oscillator, we can set 



2e^Wc 



-exp 



■0j — CjCXp 



2^2 



1 

>(0)(g+(n-l)77) 



(49) 



(50) 



where Cj (j = 1,2,3) are constants and 



2t:H 



(40) 



(41) 



Near the critical point where H = Hc2, we can neglect 
the terms since ipj £^re small. This yields the 

linearized equation, 

ai+Ki/^'^ -712 -713 \ / '^i \ 

-721 "2 + A'2/C^ -723 C2 = 0. 

-731 -732 ^3 + ^3/C^ ) \Cz ) 

(42) 



Hence, both the upper critical field Hc2 and the critical 
temperature increase as the number of bands is in- 
creased. The increase of the upper critical field, that is, 
the reduction of the coherence length ^, due to the multi- 
band effect should be noticed. In fact, the high iJc2 has 
been reported for the Fe pnictide SmFeAsOo.ssFo.isfs^. 

We can show that the three-band system is governed by 
a single Ginzburg-Landau parameter k* near i?c2- Thus 
the conventional three-band model does not yields a type- 
1.5 superconductor (37^. We define 



/ (m*)2c2/3 
2TTfi^{e*Y'' 



(51) 
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where e* = 2e, m* is the effective mass and /3 is the The free energy density is 
average of the coefficients (3j {j = 1, 2, 3). Then we can 



show, following the standard procedure[38|, |3£ 
free energy is given as 



F 
V 



1 r 
8^ 



1 + /3a(2k*2_ 1) 



that the 



(52) 



where (B^) is the integral oiBz over the space and Pa is a 
constant. This is the same formula as for the single-band 
superconductor if we substitute k to k* . The magnetiza- 
tion is similarly given by 



47rM = 



Hc2 — H 



(53) 



The multi-band effect that is contained through the pa- 
rameter K* . Hence, if there are no some reasons concern- 
ing the symmetry of crystals and superconductivity, the 
superconductor is regarded as the type 1 or type 2. 



/ = iap^ + - 2/9^ [712 COs((pi) + 723 C0S(V52) 

+ 731 COs((pi + </?2)] 



(Pn 3 



47r o 1 o 

+ ^Kp^[{V^i)^ + {V^2f + + ^2))% (58) 

where a = (1/3) and /? = (1/3) . 

The stationary conditions with respect to the fields A, 
p and (pj lead to 

c 27r / 1 27r \ 

j = rotH = 6Kc—p^ V0 + — A , (59) 
in 00 V 3 00 / 



V. CHIRALITY, KINKS AND DOUBLE 
SINE-GORDON EQUATION 

In Section III, we have shown the model that has a 
solution with time-reversal symmetry breaking. This 
section is devoted to an investigation of this solution 
in three-band superconductors by using the Ginzburg- 
Landau theory. We now consider the phase dynamics of 
the order parameters. 

A. Phase Variables 

We write the order parameters as 

^,=p,e'«^ (54) 

where pj = [tpjl is a real quantity. For simplicity, we 
assume that the coefficients of the Josephson terms are 
real: jij = 7*, — jji. The free energy density is denoted 
as /, that is, the free energy is given by the integral of / 
over the space. / is written as 

f = "j/'j ^ ^ 51 f^:>P-i ^ 2712/31/32 cos(e'i - 6*2) 
j j 

- 2723/O2P3 cos(6'2 - 6*3) - 2731P3P1 C0S(6'3 - 6*1) 

(55) 

We focus on the role of phases of the order parameters; 
we assume that 

Pj = P> I""] = (56) 
and define new phase variables 

= 01+02+^3, ^1^01-62, ^2^02-03. (57) 



6ap + 6(3p^ - 4p[7i2 cos((/7i) -I- 723 cos((p2) 

47r2 2 
+ 731 cos(<^3i)] - QKV'p + 6K—^H^x^p + -Kp{W(j)f 

00 3 

+ ^i^p[(V^l)2 + (V^2)' + (V(^2+^2))']=0, (60) 

712 sin V51-I-731 sin((/?i-f (^2)-^^[V^(/7i+V^(<^i+(/32)] = 0, 

(61) 

723sin(/3i+73isin(<y9i+(^2)-^^[V^V2+V^(<^i+</32)] = 0. 

(62) 

B. Phase Potential and Chiral States 

In this section we examine the ground state of the sys- 
tem with the potential 

V = -27i2PlP2COs((^i) - 2723P2P3C0S(((52) 

- 2731P3P1 cos((/3i + (^2)- (63) 

If we set ipz = 0'i — 01, we have ^pi + (p2 + — 
(mod27r). The minimum of this potential is dependent 
on the signs of the coefficients jijPiPj of the Josephson 
terms. We define Ti — — 2712P1P2, r2 = — 2723P2P3, and 
= — 2731P3P1. The potential is written as 

V ^Ti cos((pi) -I- F2 cos((^2) + F3 cos((^3)- (64) 

We assume that the absolute values |F,i| are almost equal 
in magnitude. Then there are four cases to be exam- 
ined as shown in Table I. When all the Fj are negative, 
we have the minimum at ipi = ip2 — = (Case 
I). If we change the sign of F3, this produces a frus- 
tration effect and ipi take fractional values. For exam- 
ple, when all the |Fi| are equal, we have a minimum at 
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where Kq = 2Kp'^ and the potential V is 
V{(p) = Vq (^cos(p + ^ cos{2ip)j 



FIG. 1: Contour map of V for Ti = r2 = Ts > 0. Black and 
white dots indicate minima of the potential V. Dotted line is 
the path in the valley connecting two minima. 



TABLE I: Classification of the ground state of the potential 
V of the Josephson interactions. Lpi in the cases III and IV 
are for IFil = IFal = irgl. 





r2 Tb ipi ip2 ip3 


Case I — 


- - 


Case II + 


+ — TT TT 


Case III - 


— + n/3 7r/3 47r/3 chiral state 


Case IV + 


+ + 27r/3 2n/3 2tt/3 chiral state 



(93i,(/32,<i23)=(7i'/3,7r/3,47r/3). In this state the order pa- 
rameters are complex and thus the time reversal symme- 
try is broken. The case IV also exhibits a similar state 
with fractional values of ipi. If all the \Ti\ are the same, 
the ground state is at {ipi,(p2,(p3)={2Tr/3,2TT/3,2n/3) 
(Fig.l). 



C. Double Sine-Gordon Equation and Kinks 

We consider the solution for 712 = 723- In this case 
we have a solution with ipi ~ ip2 = ^p. The variable tp 
satisfies the double sine-Gordon equation, 



KS/'^ip — 712 sin Lp — 731 sin(2(y9) = 0. 



(65) 



Let us consider the kink solution of the one- dimensional 
double sine-Gordon equation. The energy functional is 



E = 



dx. 



(66) 



(67) 



We defined Vq = —712^^ and u = 731/712. There are two 
cases to be examined: (1) 712 < and (2) 712 > 0. We 
show the classification of the double sine-Gordon model 
in Table II. 

First consider the case (1) Vq > 0. The potential V((p) 
has a minimum ai ip = ipQ = cos^^{—l/{2u)) if u > 1/2: 



For u < 1/2, we have a minimum at ip — n: 
n^) = Fo(-l + ^). 



(68) 



(69) 



In the case u > 1/2 we have a chiral state at ip = ipo 
and a kink solution that travels from one minimum to 
the other minimum. The double sine-Gordon equation 

§ = -;^(sin^ + ^.sin(2^)), (70) 

can be integrated for the boundary conditions dip/dx 
and ip — !• (p+ as X — >■ 00 and dtp/dx — )■ and (/? — )■ —ip+ 
a.s X —00. Here (p+ is a solution of cos{(p) = —l/{2u) 
in the range of — vr < ip < n. Since the equation above is 
integrated as 



dipV 2Vo ( " ,0 N l + 2w2\ 



we obtain the kink solution as 



ip{x) 



tan 



tan 



2b - 



1 - 



1 



'1 - 



1 

4m^ 



where q(x) is defined as 
q[x) 



exp 




(72) 
(73) 



This solution represents the one-kink or one-antikink so- 
lution. Since ip+j'K takes a fractional value, namely, 
93+ = ±27r/3 for u = 1, we call this the fractional-7r 
kink. The kink structure is shown in Fig. 2. 

In the case Vq > and u < 1/2, there are minima 
at iy9 = 7r(mod27r) in the potential V. Thus, we have a 
27r-kink solution in this case. The boundary conditions 
should be iy9 — TT as X — 00 and ip — >■ — tt as x — >■ —00, or 
vice versa. Since we obtain 



dip 
dx 



4(1- 



2 1' if 

'2u) cos — 



2u 



1 - 2-u 



(!)■ 



(74) 
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FIG. 2: (fi as a function of X = ^/2Vo/Kox for Vo > and 
u — 1, 2. Fractional kink structure is shown. For it = 1, 
ip{ — oo) — — 27r/3 and v?(cxd) — 2tt/3 which we call the 47r/3 
kink. 



FIG. 3: as a function of X = ^{u + l/2){2\Vo\/ Ko)x for 
Vo < and u — 1 and 10. This shows 27r-kink structure. For 
large u, the kink shows a saddle-like structure at a; = 0. This 
is because the potential V has a local minimum a.t ip = n for 
u > 1/2. 



the kink solution is given by 



1 



2 sinh^ (rx) 
cosh^(rx) — 2u 



(75) 



where r = ^yVo{l - 2u)/Ko. 

Second, let us consider the case (2) Vq < 0. There are 
minima at ip — (mod27r) for u > —1/2, and thus the 
2tt kink solution exists. We obtain 

/ 2 sinh^(sa;) 



2u 



1 



(76) 



Vcosh (sx) 

where s = y'lVoKl + 2u)/Kq. The kinks in this case 
are presented in Fig. 3. For large u, the kink shows a 
characteristic at a; = because the potential has a local 
minimum at (p = for u > 1/2. We have a possibility 
to find some specific features in the excited state due to 
this anomaly. For u < — 1/2 we have a fractional- tt kink 
that is given by 

'1 + 2|u|t(x)' 



(p{x) — tan 



1 



tan 



l + 2\u\/t{x) 



1 



where 




i(a:)=exp J2^H 1--^ 



(77) 



(78) 



For u = —I, this chiral solution satisfies the boundary 
condition that ip — >■ —tt/3 as x ^ —oo and ip — > 7r/3 as 
X — cx) as shown in Fig. 4. 



D. Bogomol'nyi-type Inequality 




X 



FIG. 4: ip as a. function of X = y^2\Vo\/Kox for Vb < and 
u = —1 and —2. This shows 27r-kink structure. For u — — 1, 
V5(— oo) — — 7r/3 and ifiico) = tt/3. 



and u < —1/2. Here we call the double sine- Gordon 
model with this condition the chiral sine-Gordon model. 
The energy functional for the one-dimensional system 
-L/2 <x < L/2 is 



E = 



L/2 



-L/2 



1 fdip 
2^"[d^ 



dx. 



(79) 



We assume that L is sufficiently large. If ip satisfies the 
stationary condition (5 i?/(5(^ = 0, we obtain Kod^cp/dx^ — 
dV{ip)/dip. From this we have 



The Bogomol'nyi inequality holds for the double sine- 
Gordon model with Vq > and u > l/2orVb < 



(80) 



8 



where Co = V{(po) with the condition that dip/dx 
as (/?—>■ iy9o • Then the energy is 



J {p _ 



(81) 



where = f{L/2) and = i/2), and both satisfy 
d^p/dx — > as — > For Vo > and u > 1/2, we 
obtain the energy of fractional- tt kink state as 



Ef-ktnk — 2\/2KqVqU 

+ CoL, 



1 1 
1 ^ -r^ + 71" cos 



1 

2u 



(82) 



where Co = V{ipo) = Viip+) = Viip^) = -^0(1/(4^) + 
u/2). This coincides with the energy obtained by substi- 
tuting the kink solution directly to the energy functional. 
The energy of the 27r-kink for Vb < and w > is 



= CiL + A^/ KoVa\VlT2 



1 



2V2m 



log 



2u 



1 + 2u 



/ 1 + 



2u 



1 + 2u 



(83) 



where Ci = -\Vo\{l + u / 2.) . 

Now we derive an inequahty for the energy. Let us 
consider the case Vb > and u > 1/2. By using the 
inequahty + b'^ > 2\ab\ for real a and b, we obtain 



CoL 



dx 

2 



dx 



> CoL + 2\l-KoV„u 



1 



dip 



CoL, 



dx 

sm ip-^- — smp- 



cos p - 



1 

2u 



dx 



(84) 



for the one-kink solution that satisfies dp/dx > and 
cosiy9 -f l/(2w) > 0. In the case of one fractional- tt kink 
shown above, we obtain 



E > V^KoVoi 



2W1 



1 



1 

— COS 

u 



1 

2u 



CoL. 
(85) 



where we adopt that < cos~^ {~^) — The lower 
bound of the energy coincides with the energy in eq. (|82l) . 
Here, we define the conserved current 



(86) 



with the charge 



g - / J%x)dx = —[ipix = cx^)-pix = -^)], (87) 



TABLE II: Classification of the double sine-Gordon model 
with the potential V{p) = Vo(cos((y3) + {u/2) cos(2(/9)). V{p) 
has minima ai p = po (mod27r). 

Vo u po kink 

Vb > u > 1/2 cos~^ (— l/(2w)) fractional-TT kink chiral 
Vb > u < 1/2 TV 27r-kink 
Vb < u > -1/2 27r-kink 

Vo < u < —1/2 cos^^ ( — l/(2u)) fractional-TT kink chiral 



where e^" is the antisymmetric symbol and x^ = t, 
x^ = X. A is the normalization constant defined hy A = 
cos~^(— l/(2u)) with the value in the range < ^ < tt. 
dfj_J^ ~ follows immediately from the antisymmetric 
symbol e^'^ and does not depend on the equation of mo- 
tion. The kink has Q = 1, and an antikink with Q = — 1 
exists that has a configuration with p{—oo) = A and 
p{oo) — ~A. If we assume that cosp + 1/(2m) > for 
the kink solution, the energy inequality is written as 



2KoVo 



A\Q\ 



E > sign(5^2_ftroVbw(sin(y9+ — avup. 
+ CoL. 



This is an inequality of Bogomol'nyi type. For the kink 
satisfying cosp -\- 1/(2m) < 0, the inequality is 



E > 



signQ •\/ 2Ko Vb u (sin iy9_|_ — siniy9_) 
CoL. 



2KoVr 



u 



A\Q\ 



(89) 



Since p{±oo) can take values p{oo) = ±yl-|- 2nj^T: and 
p{—oo) = ±A -|- 2n_7r with integers and n_, there 
exists a wide spectrum of solitons such that Q ~ ±.1 + 
(n+ — n_)7r/A or Q = (n+ — n-)-K/A. In general, the n- 
kink solution is a sum of (anti)kinks with plateaus where 
dp/dx vanishes. The total kink energy is a summation 
of contributions from each (anti)kink structure. Suppose 
that dp/dx = at a; = a;i, • • • , n— \ in an n-kink solution. 
We set xq = oo and x„ = -co. We define the charge 
Qj for each kink component: Qj — ip(xj-i) — p{xj). 
Q = J2i Qj holds. The energy bound is 



E > \/2KoVou'y^^ SjSigiiQj{smp{xj-^i) — smp{xj)) 



2KoVo 1 A 



u 2 



■kQ + CoL, 



(90) 



where Sj = 1 if cosp + l/(2u) > in the j-th (anti)kink 
and Sj = —1 if cos p + l/(2u) < 0, and we defined 



Q = J2sj\Qj 



(91) 
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(a) 



(b) 



(c) 



\ 71 


/ 




) 


cut 






471/3 - kink 



(-271/3,0,271/3) 



Phase difference P 
AG : = 271/3 < 



-> (271/3,0,-271/3) 



(0,271/3,-271/3) 



vortex 

^o/3 

(271/3,0,-271/3) 



271/3 - antikink 



FIG. 5: Half-quantum flux vortex with line singularity. The 
phase variables Q\ and Qi have line singularities, as shown in 
(a) and (b). Two half-flux vortices connected by the singu- 
larity are shown in (c). 



VI. FRACTIONAL VORTICES AND BOUND 
STATES 



FIG. 6: Kink, antikink and a fractional flux vortex for 
Fi = r2 = Fs > 0. The vortex is at the point P 
with flux <^o/3 where <^o is the flux quantum. We start 
from (ei, 62,613) = (-27r/3,0,27r/3) to reach (0, 27r/3, -27r/3) 
(modulo 27r) through the 47r/3-kink and 27r/3-antikink. </?i = 
di — 02 goes from —2^/3 to 27r/3 crossing the 47r/3-kink, and 
ipi goes from 27r/3 to 47r/3 = —2n/3 (mod 27r) through the 
27r/3-kink. 



In general, there are solutions of vortices with frac- 
tional quantum flux in multi-band superconductors. 
Kinks in the space of phase variables 9j play a central role 
for the existence of fractional flux vortices. In the two- 
band model, the half- quantum- flux vortex exists with a 
line of singularity of the phase variables 0j as shown in 
Fig. 5. Here 61 changes from to tt (or vr to 0) across 
the cut, and simultaneously 02 changes from to — tt (or 
— TT to 0). In the case of Fig. 5, a net-change of 9i is 
27r by a counterclockwise encirclement of the vortex, and 
that of 02 vanishes, due to singularities. Thus, we have 
a half-quantum flux vortex. This is an interpretation of 
half-quantum flux vortices in triplet superconductors (40j 
in terms of the phase of the order parameter. 

In three-band superconductors, the fractional-flux vor- 
tex exists in the chiral case as well as the non-chiral 
case. Since we have the fractional-7r kink in the chi- 
ral state (cases III and IV), the new types of vortices 
with fractional flux quanta exist on a domain wall of the 
kinkfsij. The kink considered in the previous section is 
a one-dimensional structure in superconductors. There 
are many types of kinks connecting two minima of the 
potential in three-band superconductors. 

Let us discuss the fractional vortices in the three-band 
model here. Suppose that two kinks, one is a kink and 
the other is an antikink, intersect at a point P as shown 



(-271/3,0,271/3) 




(0,271/3,-271/3) 



FIG. 7: Two-vortex bound state with line singularities in the 
time-reversal symmetry broken state. The phase variables 9i 
{i = 1,2,3) have singularities that are described by kinks in 
Fig. 6. The total flux is (j>Q. Topologically, the flux 2(j)o/3 
is equivalent to —cpo/S. Thus, this state corresponds to the 
meson under the duality transformation between charge and 
magnetic flux. 
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%/3 



-<j)o/3 



A • 

(-2jt/3,0,27t/3) 



(^o/3 



O 

(271/3,0-271/3) 

f 

B 

(0,271/3-271/3) 



• c 

(271/3,-271/3,0) 



A • 

(-271/3,271/3,0) 



2(^0 /3 



O 

/ 
/ 

(271/3,-271/3,0) 

B 

(271/3,0,-271/3) 



• c 

(0,-271/3,271/3) 

w 2(^0 /3 



FIG. 8: Three-vortex bound state with hne singularities in 
the time-reversal symmetry broken state. Each vortex has 
(/f>o/3 and the total flux is (po- The phase variables Oi {i = 
1,2,3) have singularities that are fractional-7r kinks. 



FIG. 9: Three-vortex bound state with line singularities in 
the time-reversal symmetry broken state. This state corre- 
sponds to the proton. 



in Fig. 6 in a two-dimensional xy-plane. If a vortex exists 
along the z axis just at the point P, the vortex should 
have a fractional flux quantum so that the phase change 
around the point P is 2tt. For Vq > and u > 1/2, 
this is shown schematically in Fig. 6. We set the phases 
of the order parameters (6*1,^2,6*3) — (— 27r/3, 0, 27r/3) in 
some region. After crossing the 4tt/3 kink, they become 
(27r/3, 0, — 27r/3) where the phase variables ipi and ip2 
change from — 27r/3 to 27r/3. If there is also a domain wall 
of an antikink that starts from the point P as in Fig. 6, 
we have the phases (0, 27r/3, — 27r/3) after we cross the 
antikink. Here, we obtain the phase difference between 
the initial and final states (see Fig. 6). In this case, the 
vortex that is located through the point P along the z 
axis should have a fractional flux quantum 00/3. Thus, 
in the chiral region of three-band superconductors, the 
existence of fractional vortices is easily concluded in this 
way. 

In the three-band model, the fractional vortex has two 
line singularities (kinks) in the phases of the gap function 
as shown in Fig. 6. From Fig. 6, we have a two-vortex 
bound state as presented in Fig. 7 in the chiral state. Two 
vortices form a 'molecule' by two kinks. This state may 
have lower energy than the vortex state with quantum 
flux 00 since the magnetic energy (5/9)0o smaller than 
(/)q of the unit flux. The energy of kinks is proportional 
to the distance R between two fractional vortices if R 
is large. Thus, the attractive interaction works between 
them if R is sufficiently large. 

Three-vortex bound states are also formulated: they 
are shown in Figs. 8, 9 and 10. The first two figures indi- 



cate bound states in the time-reversal symmetry broken 
state. The last one is for the unbroken state[41[. These 
states correspond to baryons if we regard the magnetic 
fiux as charge. 

VII. APPLICATIONS 

In this Section, we briefly consider some applications. 
In the one-band model, the symmetry of Cooper pairs is 
classified into irreducible representations and each rep- 
resentation is one dimensional in most cases except the 
two-dimensional representation for p-wave symmetry [42|- 
[i^ l- In the one-dimensional representation, the phase of 
the order parameter is not important, except the case 
where two states in two different representations acciden- 
tally are degenerate. In multi-band systems, the relative 
phase between bands becomes important and plays an 
essential role. 

A promising candidate of a multi-band superconductor 
is a superconductor-insulator-superconductor (SIS) junc- 
tion. An SIS tri-junction will be a three-band supercon- 
ductor, interacting through weak Josephson couplings, 
with three equivalent bands. We hope that SIS junc- 
tions or superlattice systems will be realized as artificial 
multi-band systems. 

The Fe pnictides in general contain several bands and 
it is probable that the order parameters in these bands 
belong to the same irreducible representation. We can 
expect high Hc2 due to the multiplicity of bands in Fe 
pnictides. If we assume that the pairing interactions are 
frustrating, that is, in the cases III and IV, or Vb > 
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A 

(0,0,0) 

%/3 . 



(^0/3 

o • ^ 

/ ^ \ (471/3,471/3,471/3) 

/ 

(271/3,-271/3,0) 



%/3 



B 

(27c/3,27r/3,27t/3) 



FIG. 10: Three-vortex bound state with hne singularities in 
the time-reversal symmetric state. Each vortex has <^o/3 and 
the total flux is (l>o- In this state, the region including the 
point O has higher energy. 



in the potential V, there is a possibility that chiral su- 
perconductivity will emerge in Fe pnictides superconduc- 
tors. Let us examine the pairing interactions that origi- 
nate from the spin susceptibility x(q). It has been shown 
that x(q) has peaks at q = (tt, 0) and (0, tt) in multi-band 
models for Fe pnictides |45| . In some Fe pnictides yfq) 
has a peak at q = (tt.tt) as well as at q = (7r,0)[lJ]. In 
this case we have a frustration between the pairing in- 
teraction due to x(j^j 0) and xi^^i "')■ This may lead to a 
chiral s-state ScMrai rather than s± and s+-|_ as an inter- 
mediate state of these two states. 



breaking and fractional flux vortices in the three-band 
model. Our discussion is based on the mechanism of 
time-reversal symmetry breaking due to the degeneracy 
of the gap equation. The other mechanism should also 
be investigated in future studies [4*3] . We derived the 
Ginzburg-Landau free energy for multi-component su- 
perconductors with more than two components. The 
coefficient of each term is expressed in terms of the in- 
verse of the matrix of pairing interactions G = [gij). We 
have shown that the system is determined by the single 
Ginzburg-Landau parameter k* near the upper critical 
field Hc2- For the frustrating pairing interactions, the 
chiral superconducting state exists with broken time re- 
versal symmetry. The double sine-Gordon model appears 
as an effective theory to describe low-excitation energy 
states for three-component superconductors. In the chi- 
ral case, this model has solutions of fractional- tt kinks 
and satisfies the inequality of Bogomol'nyi type with the 
topological charge Q. We also discussed the existence 
of fractional flux vortices and that of the multi-vortex 
bound state in multi-band superconductors. 



The kinks form domains, which we call the chiral do- 
mains, in multi-band superconductors; this is analogous 
to the domains in triplet superconductors f47| . In the 
chiral region of three-band superconductors, there exist 
several types of kink solitons. Vortices with fractional 
flux quanta exist on domain walls. It is important to 
investigate the effect of domain walls on the vortex dy- 
namics. Domain walls play an important role in strong 
pinning of vortices and will be closely related to flux- 
flow noise. The generation of chiral domains by quench- 
ing superconductors is also an attractive subject. This 
is a phenomenon caused by the Kibble mechanism in 
superconductors [H, |4^. In experiments, the number of 
domains may be controlled by quenching speed and sam- 
ple morphology that pins domain walls. 



VIII. SUMMARY AND DISCUSSION 

In this paper we have investigated multi-band super- 
conductors. We have examined time-reversal symmetry 
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